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1 We review the high-redshift X-ray selected AGN clustering, based on the 

\J^} ' XMM/2dF survey, and compare it with other recent XMM-based results. Us- 

, ing the recent Hasinger et al. (2005) and La Franca et al. (2005) luminosity 

, functions we find that the spatial clustering lengths, derived using Limber's in- 

■ version equation, are ~ 17 and 20 h~ 1 Mpc for the soft and hard band sources 
f - **) while their median redshifts are z ~ 1.2 and 0.8, respectively. The correspond- 

i ing bias factors at z = z are ~ 5.3 ± 0.9 and ~ 5.1 ± 1.1, respectively. Within 

f *^ , the framework of a flat cosmological model we find that our results support a 

"""■^ ■ model with Q m ~ 0.26, ag ~ 0.75, h ~ 0.72, w~ —0.9 (in excellent agreement 

' with the 3 year WMAP results). We also find the present day bias of X-ray 

Q-f AGNs to be b a ~ 2. 

o : 

+3 , 1. Introduction 

■ Active Galactic Nuclei (AGN) can be detected out to high redshifts and 
, thus their clustering properties can provide information on the large scale 

J^j ' structure, the underlying matter distribution and the evolution with red- 

shift of the AGN phenomenon. Optically selected AGN surveys miss large 
numbers of dusty systems and therefore provide a biased census of the AGN 
phenomenon. X-ray surveys, are least affected by dust providing an efficient 
tool for compiling uncensored AGN samples over a wide redshift range. 

Early studies of the X-ray AGN clustering, using Einstein and ROSAT 
data, produced contradictory results (eg.; 9 ; 43 ; 11 ; 131 ). Recently, there has 
been an effort to address this confusing issue and determine the cluster- 
ing properties of both soft and hard X-ray selected AGNs, based on the 
new XMM and Chandra missions (eg.; 46 ; 47 ; 7 , 8 ; 21 ; 36 ; 1730 ). Most of these 
studies find a relatively large correlation length for the high-z X-ray AGNs 
(eg.,; 7 ; 836 ). However, the recent XMM-COSMOS survey results 30 provided 
significantly smaller clustering lengths, opening again the contradiction 
realm. 
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Fig. 1. The XMM/2dF survey angular correlation function for the 2 bands considered. 

2. The XMM/2df survey angular and spatial clustering 

The XMM- Newton/ '2dF survey is a shallow (2-10 ksec per pointing) survey 
comprising of 18 XMM-Newton pointings equally split between a Northern 
and Southern Galactic region. Due to elevated particle background we anal- 
ysed a total of 13 pointings. A full description of the data reduction, source 
detection and flux estimation are presented in Georgakakis et al., 1918 

We derive the source log iV-log S after constructing sensitivity maps in 
order to estimate the area of the survey accessible to point sources above a 
given flux limit (see, 78 ). In table 1 we also provide the effective flux- limit 
of the different surveys, estimated using the corresponding area curves, by: 
SJj* ~ ^fiiSijm.i/^fii, where fij is the survey solid angle of which the 
flux limit is Su m ,i- 

We calculate the angular correlation function using the estima- 
tor: w(8) — /{Ndd/Ndr) ~ 1, of which the uncertainty is: a w — 
a/(1 + w(9))/Ndr, where N& d and Ndr are the number of data-data and 
data-random pairs, respectively, in the interval [8 — A8,8 + A8}. The nor- 
malization factor is / = 2Nr/(Nd — 1), with Nd and Nr the total num- 
ber of data and random points, respectively. For each XMM pointing we 
produce 100 Monte Carlo random catalogues having the same number of 
points as the real data accounting also for the sensitivity variations across 
the surveyed area (for details see, 78 ). We use a standard x 2 minimization 
procedure to fit the measured correlation function assuming a power-law 
form: w(8) = (So/d) 1 ^ 1 and fixing 7 to 1.8. The fitting is performed for 
angular separations in the range 40 — 1000 The resulting raw values of 
8 are corrected for the integral constraint and the amplification bias (for 
details see 43 and 8 ), although such corrections are quite small. The final re- 
sults are presented in Table 1. Now, the spatial correlation function can be 
modeled as (eg. 14 ): £(r, z) — (r/r )~ 7 x (l + z)~( 3+e \ where e parametrizes 
the clustering evolution model. For e = 7 — 3 (ic., e = —1.2 for 7 = 1.8), 
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Fig. 2. Le/i Panel: The expected ^-distributions using the Hasinger et al 23 soft-band 
luminosity function for the XMM/2dF survey (thick line), the XMM-COSMOS survey 
(blue long-dashed line) and the XMM-ELAIS-S1 survey (red short-dashed line). Right 
Panel: The corresponding redshift distributions using the LDDE hard-band luminosity 
function of La Franca et al. 28 The histogram in the left panel corresponds to the dis- 
tribution of the X-ray sources in the ROSAT Lochman Deep Field, 39 while that of the 
right panel on limited spectroscopic and photo-z data of our XMM/2dF survey. 



the clustering is constant in comoving coordinates, a model which appears 
to be appropriate for active galaxies (eg. 26 ). 

In order to invert the angular correlation function to three dimensions 
we utilize Limber's integral equation (eg. 34 ). For a spatially flat Universe, 
Limber equation can be written as: 

"( e >= 2 vi Udj) £ W„ £Md " (1) 

where the number of objects in a survey of a solid angle f2 s and within the 
shell (z, z + dz), is given by: 

^- = n sX ^{x){^E-\z) (2) 

with (j}(x) the selection function and x the proper distance related to the 
redshift through (see 34 ): 

X ( Z ) 4 |l With E ( Z ) = l°m(l + zf + Oa] 1 / 2 , (3) 

H J E(u) 

Since we do not have complete redshift information for our sources we 
estimate dN/dz using the X-ray source luminosity function and folding in 
the area curve, via the relation: <j)(x) — , z -> &(L X , z)dL, where $>(L X , z) 
is the luminosity dependent density evolution luminosity (LDDE) function. 
In the present analysis we use that of Hasinger et al. 23 for the soft-band 
and that of La Franca et al. 28 for the hard-band. 

In Fig. 2 we present the expected redshift distributions of the soft and 
hard X-ray sources for all three recent XMM based surveys together with 
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Table 1. Correlation function analysis results with their lrj uncertainties for the dif- 
ferent surveys and X-ray bands. The correlation lengths assume a slope 7 = 1.8. 



Band (koV) 


Survey 










M") 


r (h- 1 Mpc) 


z 


0.5-2 


XMM/2dF 


432 


1.3 


xlO" 


14 


10.4 ± 3 


17 ± 2.2 


1.25 


0.5-2 


XMM-COSMOS 


1037 


1.9 


xlO" 


15 


1.9 ± 0.3 


9.3 ± 0.7 


1.33 


0.5-2 


XMM-ELAIS 


395 


4.5 


xlO" 


15 


5.2 ± 3.8 


13.8 ± 4.5 


1.32 


2-8 


XMM/2dF 


177 


4.6 


xlO" 


14 


22.2 ± 9 


20 ± 3.7 


0.80 


4.5-10 


XMM-COSMOS 


151 


1.6 


xlO" 


14 


6 ± 2 


12 ± 1.8 


0.89 


2-10 


XMM-ELAIS 


205 


1.8 


xlO" 


14 


12.8 ± 7.8 


17 ± 4.6 


0.88 



the histogram of some limited spectroscopic and photo-z data (see caption 
for details) . It is evident that all three surveys trace similar depths although 
the XMM-COSMO and XMM-ELAIS surveys have a slightly larger contri- 
bution from zj> 2 in the soft-band. The predicted median redshift of the 
three XMM surveys for the soft and hard sources are shown in table 1. 

Then the inversion of eq.(2), using the LDDE luminosity function, e = 
— 1.2, the concordance cosmological model, and integrating out to z = 4.5, 
provides a correlation length of r ~ 17 ± 2.2 h^ 1 Mpc and ~ 20 ± 3.7 h^ 1 
Mpc, for the soft and hard bands, respectively (the slight differences with 
respect to the Basilakos ct al., 78 results are due to the different luminosity 
functions used and a better treatment of the errors) . These results are in 
good agreement with the XMM-ELAIS results 36 but significantly larger 
than the XMM-COSMOS 30 results (see table 1). 

If we assume that all three XMM based surveys provide correct re- 
sults, and that their differences are due to the cosmic variance, we can 
volume weight them to obtain an average estimate of the clustering length 
of the XMM point sources. The relative volume weights for the XMM- 
COSMOS, XMM-ELAIS and XMM/2dF surveys are 1:0.9687:0.8527 and 
1:0.9947:0.701 for the soft and hard bands, respectively (we have not cor- 
rected for the slight differences of the hard-bands used in the three sur- 
veys). Doing so we obtain: (r ) — 13.2 and 16.0 h^ 1 Mpc, for the soft and 
hard bands, respectively These large correlation lengths are comparable 
to those of Extremely Red Objects, 20 of luminous radio sources ( ; 37 ; 3238 ) 
and of bright distant red galaxies. 15 They are however, significantly larger 
than those derived from optical AGN surveys: r ~ 5.4 — 8.6 h^ 1 Mpc 

f„. .12.27.13.22.3544\ 
\eg., , , , , ). 
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3. Cosmological Constraints 

It is well known 24 that according to linear biasing the correlation function 
of the mass-tracer (£ bj) and dark-matter one ((dm), are related by: 



where b(z) is the bias evolution function. In this study we use the bias 
model of Basilakos & Plionis, 34 We quantify the underlying matter distri- 
bution clustering by presenting the spatial correlation function of the mass 
£dm(P) z) as the Fourier transform of the spatial power spectrum P(k): 



where k is the comoving wavenumber and e = —1.2, according to the con- 
stant in comoving coordinates clustering evolution model. As for the power 
spectrum, we consider that of CDM models, where P(k) = P k n T 2 (k) 
with scale- invariant (n = 1) primeval inflationary fluctuations. In partic- 
ular, we use the transfer function parameterization as in, 2 with the cor- 
rections given approximately by Sugiyama. 41 Note that we also use the 
non-linear corrections introduced by Peacock & Dodds. 33 We have cho- 
sen to use either the standard normalization given by: as — 0.5fi m 7 with 
7 ~ 0.21 — 0.22w + 0.33il m , 45 or to leave eg a free parameter. 

Firstly, using equations (1), (4) and (5) and evaluating at z = z we 
derive the bias factor of our sources at the corresponding redshift. We find 
b(z) ~ 5.1 ± 1.1 and ~ 5.3 ± 0.9 for the hard and soft bands respectively, 
which are significantly smaller than the values attributed to our work in 
Miyaji et al. 30 

Next, we utilize a x 2 likelihood procedure to compare the measured 
XMM soft source angular correlation function with the prediction of dif- 
ferent spatially flat cosmological models (see, 56 ). In particular, we define 
the likelihood estimator as: £ AGN (c) oc cxp[— Xagn( c )/2] with: Xagn( c ) = 
zC"=i [ w th(6i, c) — w b s (9i)/<Ji] 2 , where c is a vector containing the cosmo- 
logical parameters that we want to fit and &i the observed angular corre- 
lation function uncertainty. We assume a flat (O to t = 1) cosmology with 
primordial adiabatic fluctuations and baryonic density of fl^h 2 ~ 0.022 
(eg.; 2540 ). In this case the corresponding vector is c = (f2 m , w, cr 8 , h, b ) 
and we densely sample the various parameters. Note that in order to inves- 
tigate possible equations of state, we allow the parameter w to take values 
below -1, corresponding to the so called phantom cosmologies (eg. 10 ). 

The resulting best fit parameters for e = —1.2 and the Miyaji et al 29 
luminosity function are presented in Table 2. In the first two rows we present 



€obj(r,z) = b 2 (z)£ DM (r,z) 



(4) 
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Table 2. Cosmological parameters from the likelihood analysis. Errors of the 
fitted parameters represent la uncertainties. Note that the fitted parameters cor- 
respond to results marginalized over the parameters that do not have errorbars, 
for which we use the values indicated. 



Data 




C8 


w 


h 


bo 


XMM 


qi+ 016 


0.93 


w — —1 


n 72+ 02 


2 S0+ - 70 
z - ou -0.20 


XMM/SNIa 


0.28 ± 0.02 


0.95 


-1 05+ 010 
i.uo_ 20 


0.72 


2.30 


XMM 


0.28 ± 0.03 


0.75 ± 0.03 


w — —1 


0.72 


2 -0io.25 


XMM/SNIa 


0.26 ± 0.04 


0.75 


n q+ 010 

""'-0.05 


0.72 


2.0 



results based on the traditional 45 as normalization. Note that our estimate 
of the Hubble parameter h is in very good agreement with those derived by 
the HST key project. 16 In the last two rows of Table 2 we leave as free but 
fix the Hubble constant to h = 0.72. In this case our fit provides a a s value 
which is in excellent agreement with the recent 3-years WMAP results 40 
Therefore, allowing values w < — 1 we can derive a (f2 m , rr 8 ) relation, a good 
fit of which is provided by : 

a 8 = 0.34(±0.01) ft m 7(^,«0 ( 6 ) 

with7(f2 m ,w) = 0.22(±0.04)-0.40(±0.05)w-0.052(±0.040)n m . Note that 
w is degenerate, within the la uncertainty, with respect to fl m . Therefore, 
in order to put further constraints on w we additionally use a sample of 
172 supernovae SNIa of Tonry et al. 42 

The joined likelihood analysis is performed by marginalizing 
the X-ray clustering results over a 8 , h and b - The vector c 
now becomes: c = (£l m ,w). The SNIa likelihood function can 
be written as: £ SNIa (c) oc exp[-Xg NIa (c)/2], with: x^ NIa (c) = 

zCili [log-^L 1 ( z i 7 c ) — log!?L bs ( 2; i)/ cr i] ' w hcrc D-l(z) is the dimensionless 
luminosity distance, D^(z) = H Q (1 + z)x(z) and Zi is the observed red- 
shift. The joint likelihood function peaks at: f2 m = 0.26 ± 0.04 with 
w = —0.90^0 ^5- Using eq.(6) we find that the normalization of the power 
spectrum that corresponds to these cosmological parameters is as — 0.73, 
in excellent agreement with the recent 3-year WMAP results. 40 
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